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Host galaxies and black hole masses of low and high luminosity 
radio loud active nuclei. 
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ABSTRACT 

We investigate the host galaxy luminosities of BL Lac Objects (BLLs) and Radio Loud 
Quasars (RLQs) at z<0.5 imaged with the Hubble Space Telescope (HST). From a homo- 
geneous treatment of the data we construct the host galaxy luminosity functions (HGLFs) 
and find that RLQ hosts are ~ 0.5 mag brighter than those of BLL: < Mr >rlq = -24.0, 
< Mr >bll= -23.5. For both classes the HGLFs exhibit a remarkably different distri- 
bution with respect to that of normal (inactive) ellipticals, with clear preference for more 
luminous galaxies to show nuclear activity. We make use of the black hole mass - bulge lumi- 
nosity (Mbh -L&uige) relation, derived for nearby inactive ellipticals, to estimate the central 
black hole mass in our sample of radio loud active galaxies. In spite of a ~ 2 order of mag- 
nitude difference of intrinsic nuclear luminosity BLL and RLQ have BH of similar mass 
(< Mbh/Mq >bll=5.6x 10 8 , < Mbh/Mq >r L q= LOx 10 9 ). This implies that the two 
types of objects are radiating at very different rates with respect to their Eddington luminosity. 
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1 INTRODUCTION 

There is a general consensus about the existence of supermassive 
black holes (SBH) at the center of normal galaxies as well as in the 
nuclei of active galaxies and quasars (see e.g. the recent review of 
Ferrarese 2002). A large body of data, in particular based on high 
resolution HST observations, is now available to strongly support 
the presence of such massive BH. 

SBHs play an important role in the formation and evolution 
of massive galaxies and are also a key component for the develop- 
ment of the nuclear activity. In spite of this apparently ubiquitous 
presence of SBH in galaxies our understanding on how the galax- 
ies and their central BHs are linked in the process of formation of 
the structures remains unclear but several attempts of explanation 
have been proposed (e.g. Silk & Rees 1998; Haehnelt & Kauffmann 
2000; Adams et al 2001; Burkert & Silk 2001; Merritt & Ferrarese 
2001a; Balberg & Shapiro 2002). 

From the observational point of view it was shown that the 
BH mass is correlated with the properties of the bulge compo- 
nent of the host galaxy, which is translated into the relationships 
between Mbh and the bulge luminosity (Kormendy & Richstone 
1995; Magorrian et al. 1998; Richstone et al. 1998; Kormendy & 
Gebhardt 2001) and between Mbh and the velocity dispersion a 
of the host galaxy (Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Merritt & Ferrarese 2001b). Any theory of SBH and galaxy for- 
mation must therefore take into account and explain such observed 
empirical relations (e.g. Silk & Rees 1998, Haehnelt & Kauffmann 



2000; Ciotti and van Albada 2001). On the other hand, although 
these relations have a significant scatter ( ~0.4 in Log Mbh), they 
offer a new tool for evaluation of BH masses in various types of 
AGN if a reliable measurement of the host galaxy luminosity or 
of a is done. While for AGN with strong emission lines (as QSO 
and Seyfert galaxies) the standard methods (e.g. reverberation map- 
ping) under virial assumptions of the emitting regions can be used 
to derive Mbh ( see e.g. Wandel et al. 1999 and Kaspi et al. 2000), 
the above relations may be the only way to estimate Mbh for ac- 
tive galaxies that lack of emission lines or that are too far away (as 
BL Lac objects and many nearby radiogalaxies) to resolve the re- 
gion of influence of the BH with present-day instrumentation. The 
two different approaches lead to consistent estimates of BH masses 
within the assumed uncertainties of the two methods (McLure & 
Dunlop 2002). 

Given the difficulty to obtain a from spectroscopy of the 
galaxies hosting active nuclei, only for few AGN it was possi- 
ble to use a to evaluate Mbh (Ferrarese et al. 2001, Barth et al. 
2002,2003; Falomo et al. 2002,2003). On the contrary the galaxy 
luminosity is much easier to measure for active galaxies and can 
therefore be used to determine Mbh for larger data set. 

In this paper we use the Mbh - Lbuige relation to investi- 
gate and compare the BH mass distribution of a sample of low 
and high luminosity radio-loud AGNs (BL Lacs and RLQs respec- 
tively). Both classes are found to reside in massive giant ellipticals 
(Urry et al. 2000; Dunlop et al. 2001), which makes them rather 
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homogeneous for such kind of analysis. To ensure uniformity of 
the results we have considered only objects at z<0.5 and that have 
been imaged by HST. This allows us also to better constrain their 
host properties. In section 2 we describe our samples of BLLs and 
RLQs and compare their host galaxy luminosity functions. In sec- 
tion 3 we discuss the MBH-Lbuige relation and derive the central 
black hole mass for each object. Finally in section 4 we discuss 
our findings comparing with recent results on radiogalaxies. In our 
analysis H o =50 Km s _1 Mpc -1 and S1 =0 were used. 



2 LUMINOSITY OF THE HOST GALAXIES 

We have collected host galaxy data for BL Lacs and RLQs at z<0.5 
imaged by HST with the WFPC2 and have constructed a homoge- 
neous dataset of the host galaxies luminosities. This yields a sample 
57 BL Lacs and 1 8 RLQs that represent, respectively, low and high 
luminosity radio loud active galaxies. 

Since most of the observations were obtained in the F702W 
and F675W filters we converted all HST magnitudes into R Cousins 
band (Holtzman et al. 1995). In the few cases where the filters 
F555W, F606W were used we applied a color correction V-R = 
0.61 for the elliptical host galaxies (Fukugita et al. 1995). Absolute 
magnitudes have been k-corrected following Poggianti (1997) pre- 
scriptions and corrected for galactic reddening using the Bell Lab 
Survey of neutral hydrogen N_s (Stark et al. 1992) with the con- 
version logN w /E(B-V)=21.83 cirT 2 mag -1 (Shull & van Steen- 
berg 1985) assuming a total-to-selective extinction A.r=2.3E(B-V) 
(Cardelli, Clayton & Mathis 1989). Since the objects are distributed 
over a significant redshift interval we have also applied a correction 
to set the host galaxy luminosity to present epoch assuming a pas- 
sive stellar evolution for massive ellipticals (Bressan et al. 1994). 
This correction (Am ~ -0.2 ) allows us to properly use the Mbh- 
Lbuige relations which refers to local galaxies. In the following M_r 
represents the host galaxy absolute magnitude including all correc- 
tion terms specified above. 



2.1 The BL Lac Objects sample 

The HST snapshot image survey of BL Lacs (Urry et al. 2000, 
Scarpa et al. 2000) has provided a homogeneous set of 1 10 short ex- 
posure high resolution images through the F702W filter. From this 
we have extracted all resolved objects at z<0.5, yielding 57 sources 
with redshift between 0.027 and 0.495 « z >=0.20±0.1 1). The 
host galaxy morphology of these objects is always well described 
by an elliptical model (Scarpa et al. 2000). The absolute Mr mag- 
nitude for each object is reported in Table 1. The host galaxy av- 
erage luminosity is < Mr >=-23.49±0.5, roughly one magni- 
tude brighter than the characteristic galaxy magnitude M R = -22.75 
(Metcalfe et al. 1998). According to the shape of their spectral en- 
ergy distribution (SED), BL Lacs are broadly distinguished into 
two types (see Padovani & Giommi 1995) : those the SED of which 
peaks at near-infrared/optical and the 7-ray MeV regions (low fre- 
quency peaked BL Lacs or LBL), and those that have SED peaking 
in the UV/X-ray and the 7-ray TeV energies (called high frequency 
peaked BL Lacs or HBL). As shown by Urry et al. 2000 the host 
galaxy properties of HBL and LBL Lacs are indistinguishable, and 
therefore the two subclasses will not be separated for this analysis. 



Table 1. Host Galaxies Properties of BL Lac Objects: (a) name of the 
source, (b) redshift, (c) absolute R host galaxy magnitude, (d) black hole 
mass in units of Mq in logarithmic scale. For reference see Urry et al. 
2000, Scarpa et al. 2000. 
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Figure 1. Redshift Distribution for RLQ (a), BLL (b) and BLL matched 
samples (see text) (c). 



Table 2. Average properties of three subsample of RLQs: (a) The sample 
(BK=Bahcall et al. 1997 and Kirhakos et al. 1999; Bo=Boyce et al. 1998; D: 
Dunlop et al. 2001); (b) Number of objects; (c) average redshift; (d) average 
absolute host galaxy magnitude. 

Sample N < z > < Mr > 
(a) (b) (c) (d) 



BK 8 0.26 -23.92±0.61 

Bo 5 0.30 -24.23±0.27 

D 10 0.22 -24.02±0.29 



2.2 The RLQs sample 

Since there is not a homogeneous and large set of HST observa- 
tions for RLQs, we have constructed a sample of 18 RLQ from 
the merging of three different subsets (BK: Bahcall et al. 1997 
and Kirhakos et al. 1999; Bo: Boyce et al. 1998; D: Dunlop et 
al. 2003). Bahcall et al. 1997 and Kirhakos et al. 1999 studied 8 
RLQ in the F606W and F555W filters and in the redshift range 
0.158< z <0.367; Boyce et al. 1998 reported the analysis for 5 
sources with 0.223< z <0.389 in the F702W filter. The largest 
subsample was investigated by Dunlop et al. (2001), who report 
host galaxy properties for 10 radio-loud quasars with 0.1< z <0.25 
observed in the F675W filter. As in the case of BLL an elliptical 
model is always a good representation for the host galaxies. The 
average properties of the three subsamples are reported in Table 
2. Our evaluations of Mr are consistent with absolute values re- 
ported by the quoted authors when galactic extinction, filter cor- 
rection and evolution correction are taken into account. Since these 
subsets have statistically indistinguishable host luminosity distri- 
butions we have merged these subsamples to construct a represen- 
tative sample of RLQ (see also Treves et al. 2002), taking average 
values for objects observed twice. The combined dataset consists 
therefore of 18 objects with redshift in the range 0. 158< z <0.389, 
< z >=0.26±0.07 and < Mr >= -24.04 ± 0.4. In Table 3 we 
give for each source the redshift z, the assumed galactic extinction 
in the R-band Ar, the host galaxy apparent magnitude R and the 
absolute magnitude Mr. In figure 1 we compare the redshift distri- 
butions for the RLQ and BLL samples. 



Table 3. Host Galaxies Properties of RLQs: (a) name of the source, (b) 
redshift z, (c) adopted extinction in the R band Ar, (d) apparent host galaxy 
R magnitude, (e) absolute R host galaxy magnitude, (f) black hole mass in 
units of M0 in logarithmic scale and (g) reference, where BK=Bahcall et 
al. 1997 and Kirhakos et al. 1999, Bo= Boyce et al. 1998, D=Dunlop et al. 
2003. 
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2.3 Comparison of the host luminosity of BLL and RLQ 

Both BLL and RLQ in our sample have been mostly discovered as 
counterparts of radio and/or X-ray sources. Therefore the objects 
considered here were selected on the basis of the nuclear proper- 
ties, and because there is not a significant correlation between the 
nuclear and host galaxy luminosity (Urry et al. 2000; Percival et al. 
2001; Dunlop et al 2003 ), we can consider the distribution of the 
host galaxy luminosity is unbiased by selection effects. Moreover 
the homogeneous treatment of the data attests a reliable comparison 
of host luminosity between the two classes (BLLs and RLQs). The 
main limitation of this comparison remains however the exiguity of 
the RLQ sample. 

We find that the average absolute magnitude of RLQ is about 
0.5 magnitude brighter than that of BLL. The difference is illus- 
trated in figure 2 where we compare the cumulative absolute mag- 
nitude distributions of the hosts for the two samples (a KS test indi- 
cates that they are statistically different at the > 99% level). Since 
the two samples span slightly different redshift range we checked 
that this does not affect our result. If we consider a subsample of 
BLL with redshift distribution matched with that of RLQ (see fig- 
ure 1 (c)) we find < Mr >bll (matched)= -23.54±0.47 thus con- 
firming our finding. Given the homogeneity of data analysis and the 
procedure for the selection of the objects we believe that this dif- 
ference is not biased. We note, however, that a larger number of 
objects (in particular of RLQ) is required to confirm this result on 
a firm statistical basis. 

To further compare the luminosity distributions of the host 
galaxies we constructed the host galaxy luminosity function 
(HGLF) for the two subsets of objects. To set the normalization 
of the HGLFs we simply assume the space density of both class of 
objects as derived from studies of complete samples. For BL Lacs 
we use the value <J?o= 10 Mpc~' mag -1 of the FR I radiogalax- 
ies luminosity function at Mh=-22.8 given by Padovani and Urry 
(1991) under the assumption that FR I radiogalaxies are the parent 
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Figure 2. Comparison of the cumulative host galaxy absolute magnitude (R 
band) distributions of RLQ (dotted line) and BL Lacs (solid line). 

population of BL Lacs (e.g. Urry and Padovani 1995). For RLQ, 
we took the value of the LF of close-by radio quiet QSO (Koehler 
et al 1997, Grazian et al 2000) at Mb= -25.1, which corresponds 
to the average value of the nuclear magnitude for our sources, and 
scaled it by a factor 10 to account for the ratio between RQQ and 
RLQ (e.g. Moderski et al 1998). This yields $ =2.3 x 10 -9 Mpc" 3 
mag -1 . 

In figure 3 we show the HGLF of BLL and RLQ compared 
with that of inactive ellipticals (Metcalfe et al. 1998). To quantify 
the differences in shape of the HGLF we fitted the luminosity dis- 
tributions of the host galaxies with a modified Schechter function 
$=K x $s x (L/L*)' 3 , where <J>s is the Schechter function for el- 
liptical galaxies (Metcalfe et al. 1998): $ s =$* x (L/LT x exp(- 
L/L*), assuming $* = 8.5 x 10~ 2 Mpc -3 , a=-l.2 and L*=2.25 
x 10 44 erg s -1 (Metcalfe et al. 1998). The best fit to HGLF was 
estimated minimizing \ 2 f° r the function $. We find /3 =2.7 ±0.2 
for BLLs, p =3.6±0.3 for RLQ. The shapes of the two HGLF are 
somewhat different, but only at the 2-a level. 

This suggests that a given elliptical has a probability of having 
a radio-loud active nucleus depending on the galaxy luminosity. 
Moreover one can argue that the steepness of this behavior depends 
on the intrinsic luminosity of the nucleus as hinted by the different 
value of p for BLL and RLQ. 

It turns out therefore that both types of radio loud active galax- 
ies exhibit a remarkable different distributions with respect to nor- 
mal ellipticals, with clear preference for more luminous (and mas- 
sive) galaxies to show nuclear activity. This behavior disagrees with 
that found by Wisotzki et al. (2001) for the host galaxies of radio- 
quiet QSO. The shape of HGLF of the latter objects in fact appears 
to be consistent with that of ordinary inactive early type galaxies. 



3 MASS OF THE CENTRAL BLACK HOLE 

Basing on dynamical studies of nearby early type galaxies it was 
shown that there is a linear relation between the luminosity of the 
spheroidal component of a galaxy (Lbuige) and the mass Mbh of 
the central black hole (e.g. Kormendy & Gebhardt 2001 and refer- 
ences therein). This correlation has a scatter of ~ 0.4 in Log(MBH ) 
that can be ascribed mainly to the errors of measurements of Mbh 
and to the uncertainties to disentangle the bulge from the disc com- 
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Figure 3. The HGLF of RLQs and BL Lacs {filled points) compared with 
the fit (solid lines ) with a modified Schechter function (see text). A slightly 
different value of /3 is found for BLL (J3 =2.7) and RLQ (/3 =3.6). The 
dotted line is the fit to RLQ data normalized to the BLL data. The dashed 
curve represents the luminosity function of elliptical galaxies of Metcalfe 
etal. (1998). 

ponent of the galaxies. Nevertheless it can be used to estimate Mbh 
for our objects provided that consistent and reliable host galaxy lu- 
minosities are used. 

To avoid systematic effects it is important that the adopted ab- 
solute magnitude of the galaxy be homogeneous (in terms of spec- 
tral band, adopted cosmology, extinction correction, filter, etc) with 
that used to derive the Mbh -Lbuige relation. To satisfy this require- 
ment we used the relationship between Mbh and M_r derived by 
Bettoni et al. (2003) for 20 inactive ellipticals assuming the same 
calibrations of Mr we have adopted here. 

From the analysis of the inactive sample of ellipticals Bettoni 
et al. give 

log{MsH I Mq ) = -0.50 x M R - 3.00 (1) 

that is used to derive Mbh from absolute (total) magnitude 
Mr of ellipticals (Ho=50 Km s -1 Mpc -1 ). This relation has rms 
scatter of 0.38 in log(Msn) and it is similar to that derived by 
McLure & Dunlop (2002). 

Since the host galaxies in the samples considered here are all 
bona-fide ellipticals we used relation (1) to derive Mbh for the two 
samples of radio loud AGN (BL Lacs and RLQs) and report the 
value for each object in table 1 and 2, respectively. The uncertainty 
on the estimated black hole mass is dominated by the scatter of 
relation (1) while the uncertainties on the host galaxy magnitude 
are usually smaller. 

The distributions of Mbh for the two samples are shown 
in Figure 4. The two classes exhibit an average difference by 
a factor ~ 2 in Mbh as a consequence of the different aver- 
age host luminosity. We find that the average values of Mbh are 
< log (Al b h / Mq) >= 8.75±0.25 and < log(M B H/M Q ) >= 
9.02±0.20 respectively for BLLs and RLQs. 

A complementary method to derive BH masses is based on the 
relation between the BH masses and a of the host galaxy (Gebhardt 
et al. 2000 and Ferrarese et al. 2000). This has been applied to a 
small number of nearby BL Lac objects (Falomo et al. 2002,2003; 
Barth et al. 2002;2003). As shown by Falomo et al 2003 there is 
a good agreement in the results obtained with the two techniques. 
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Figure 4. Comparison of the Black Hole Mass (derived from the host galaxy 
luminosity using the relation between Mjg and bulge luminosity of nearby 
ellipticals) distributions of RLQs and BLLs 



Thus no estimates of BH mass of RLQ are available from a because 
of the lack of measurements. 

Wu, Liu and Zhang (2002) and Woo and Urry (2002) have de- 
rived estimates of the BH mass of BL Lacs using the Mbh - u re- 
lation where the latter quantity was inferred from measurements of 
the effective surface brightness and the effective radius of the host 
galaxy and by assuming these are linked to a through the Funda- 
mental Plane relationship. Although in principle this method could 
work we believe it is less accurate than the direct use of the Mbh 
- M(bulge) relation. In fact in addition to the uncertainty in the 
measured quantities (/i e and R e ), which are much larger than the 
total magnitude of the galaxy, one has to take into account the un- 
certainty due to the intrinsic scatter of the FP relation and that of 
the Mbh - o relation. Comparison of our BH masses (see Table 1) 
with those derived by Woo and Urry or Wu et al 2002 indicate re- 
markable differences of Mbh for many objects mainly because of 
the wrong evaluation of the velocity dispersion. In some cases (e.g. 
1807+698, 1104+384) the poor estimate of the velocity dispersion 
via the FP method has been clearly confirmed by direct measures 
of a (Barth et al; Falomo et al). 

We also note that in the Woo & Urry (2002) estimates addi- 
tional errors derive from a mistreating of the FP parameters because 
instead of the average surface brightness (< fi > e ), the isophotal 
surface brightness was used. 



4 DISCUSSION 

The analysis of HST images for low redshift BL Lacs and radio 
loud quasars has shown that for both types of active nuclei the host 
galaxies are very luminous ellipticals. On average they are ~l-2 
mag more luminous than the typical galaxy luminosity (M|j ~ - 
22.75; Metcalfe et al. 1998). After homogeneous treatment of the 
data we also found that host galaxies of RLQ are systematically 
more luminous by ~0.5 magnitudes than BL Lac hosts. Although 
this result does not seem to depend on the selection of the objects, 
a larger sample of RLQ is needed to reach a firm conclusion. We 
have shown that the distribution of the host galaxy luminosity ex- 
hibits a marked drop towards less luminous galaxies and that it is 
somewhat different for the two classes (BLL and RLQ). This in- 



dicates that such kind of nuclear activity occurs preferentially (or 
lasts longer) in massive galaxies. The different distributions of host 
luminosity for BLL and RLQ may simply reflect the very large 
range of intrinsic nuclear luminosity (about two orders of magni- 
tude, see below). High power nuclear activity like that observed in 
RLQ can occur only in the most luminous and massive galaxies and 
it is therefore a rare event. On the other hand, low power nuclear 
activity as that observed in BL Lacs (or in radio galaxies believed to 
be identical objects non affected by beaming effects) can be present 
also in galaxies with intermediate luminosities. 

On the assumption that the galaxy luminosity is correlated 
with the central BH mass, the host galaxy luminosity can be trans- 
lated into central BH masses. It turns out thus that, within a factor 
of two, BLL and RLQ have similar BH masses but their total in- 
trinsic nuclear luminosities are remarkably different. In addition 
to the higher observed nuclear/host ratio of RLQ with respect to 
BLL we have to take into account the fact that, while we consider 
RLQ basically unbeamed, for BLL a substantial beaming factor is 
present (S ~15 see Ghisellini et al. 1998, Capetti & Celotti 1999). 
The intrinsic nuclear luminosities therefore differ by about a fac- 
tor 100. This implies a dramatic difference of the Eddington ratio 
£ E = L/L E where L E = 1.25 x 10 38 x(M bh /Mq) erg s _1 (see 
also O'Dowd et al 2001, 2002; Treves et al 2002). Basing on the 
estimated total QSO luminosity of L~ 3x 1O 12 L (e.g. Elvis et 
al. 1994) and assuming BH masses of 1-5 x 10 9 Mq, we find that 
RLQ may be emitting at rates of 10% or higher than their Edding- 
ton power, while BLL are always emitting at regimes that are much 
lower than L E . 

According to the unification schemes of radio loud AGN (e.g. 
Urry and Padovani 1995) BL Lacs are radiogalaxies the jet of 
which is closely oriented toward the observer. Based on arguments 
of number density, luminosity functions and unbeamed properties 
(as the extended radio luminosity or the host galaxies) the par- 
ent population of BL Lacs is likely formed by FR I radio galax- 
ies with some contamination by FR II sources (Padovani and Urry 
1990; Wurtz et al 1996; Falomo and Kotilainen 1999; Cassaro et 
al 1999; Urry et al 2000). Under this hypothesis the BH mass of 
BL Lacs « log(M B H/M ) >= 8.75±0.25 ) and of the parent 
(unbeamed) objects must be identical. In figure 5 we compare the 
distribution of the BH mass for our sample of BL Lacs with that 
of low redshift radiogalaxies from the sample of Govoni et al 2000 
(see also Bettoni et al 2003). The two distributions are rather sim- 
ilar although the most massive BH in luminous FR I radiogalaxies 
do not appear to have counterparts in the known BL Lacs. The av- 
erage values of Mbh are < log(M B H/M e ) >= 9.04±0.30 and 
< log(M B H/M Q ) >= 8.78±0.35 for, respectively, FR I and FR 
II radio galaxies. 

Finally we wish to note that in addition to the mass the other 
parameter which characterizes a black hole and that may play a rel- 
evant role in the observed phenomenology is the BH spin. It has 
been suggested that the spin energy is responsible for the jet emis- 
sion hj and therefore for the development of the radio-emission 
(e.g Blandford 2000, Dunlop et al 2003). The spin is clearly not 
directly measurable but it could be deduced from an estimate of Lj 
and the BH mass using the Blandford and Znajek (1977) formula. 
While Lj could be obtained from the spectral energy distribution 
(e.g. Tavecchio et al 2000, 2002), the BH mass may come through 
the procedures described in this work. 
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Figure 5. Black hole mass distribution of our sample of BL Lacs (top panel) 
compared with that of FR I (middle panel) and FR II (bottom panel) radio- 
galaxies studied by Bettoni et al (2003). The open histograms in the middle 
and bottom panels refer to the distribution of the whole sample of radio- 
galaxies (see text). 
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